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We study the dynamics of a single excitation in an infinite XXZ spin chain, which is launched 
from the origin. We study the time evolution of the spread of entanglement in the spin chain and 
obtain an expression for the second order spatial moment of concurrence, about the origin, for both 
ordered and disordered chains. In this way, we show that a finite central disordered region can lead 
to sustained superballistic growth in the second order spatial moment of entanglement within the 
chain. 
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A. Introduction 

Entanglement can be viewed as a physical resource 
that has no analog in classical information theory. As 
such, entanglement plays an important role in many 
quantum information tasks, such as quantum cryptogra- 
phy, teleportation and quantum algorithms. The devel- 
opment of protocols for the distribution of entanglement 
is an important problem in quantum information pro- 
cessing. Recently several methods have been proposed 
for accomplishing the relatedproblem of quantum state 
transfer using spin chains ELIII1I1I1I1I1I1I1I3I3 
13 1X3- • An important characteristic in all of these 
models is that the transfer occurs with minimal control 
or intervention in the dynamics of the spin chain. Besides 
achieving perfect quantum state transfer, another impor- 
tant task that can be examined is that of distributing en- 
tanglement between a central node and many distributed 
nodes within a processor. Such a capability may be useful 
for the generation of multi-qubit entangled states (with 
the appropriate post-distillation). In any case, gaining 
a better understanding of how one can use spin chain 
dynamics to spread out (rather than transmit) entangle- 
ment is of fundamental interest, although this has not 
seemed to have been addressed much in the literature to 
date. 

In this letter we consider the second order spatial mo- 
ment of entanglement in a spin chain, with a single exci- 
tation, and show that a finite disordered region can lead 
to superballistic growth in the the second order spatial 
moment. The entanglement properties of various spin 
chains have been examined in a number of recent works 
CI 13 13 El 13 E3- It has been shown that, for 
a particle in a lattice with a finite disordered region, a 
finite period of superballistic growth can occur in the 
spatial variance of the wave function [2lJ. This can be 
applied to XXZ spinchains, which have nearest neighbor 
interactions similar to the tight binding approximation 
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FIG. 1: The second order spatial moment of concurrence 
M(t), for an ordered (dark grey) and a disordered (light grey) 
chain. Each chain is 40500 sites long and the disordered region 
in the second chain is 100 sites long. The J z couplings in the 
central disordered region vary randomly, uniformly between 
and 2.5. A period of super-ballistic, ~ t 3 , spatial expansion 
of the moment can be observed, for a period, which tends to 
a long time super-ballistic expansion rate of ~ t 5 ^ 2 , exceeding 
the long time expansion rate of the ordered chain, ~ t 2 . 

used for a 1-dimensional lattice. Unlike the spatial vari- 
ance, however, the growth in the second order spatial 
moment of concurrence in a chain with a finite central dis- 
ordered region is not bounded from above by the similar 
growth experienced in an ordered chain. The disordered 
region leads to sustained superballistic spatial growth of 
the concurrence and provides an efficient mechanism for 
distributing entanglement along spinchains. 

B. Concurrence 
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One measure of bipartite entanglement is concurrence 
[2^ . which is a monotonic function of entanglement of 
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formation. For a mixed state, the concurrence between 
two qubits, i and j, is defined to be = max{Ai — 
A2 — A3 — A4,0}, where A„ is the square root of the n th 
eigenvalue of pp in descending order. Here p = (a y (S> 
a y )p*(a y £g> <J y ), where p* is the complex conjugate of p. 

For a single pure excitation of a N— spin 1/2 chain, the 
wave function is 



= ai \l) + ■ ■ ■ + a N \N) 



(1) 



where \k) is one of N basis states of the A^— spin chain 
where the k th site is in the excited state and all other sites 
are in the ground state. Taking any two spins within this 
chain i, j, and tracing over the rest yields the following 
density matrix: 



where p, = 1 
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The eigenvalues of PijPij are 



j 1 2 , 0,0,0}, so 

A„ = {2| ai 1 1 |,0, 0,0}. Inserting these expressions for 
A„ into the previous equation for concurrence gives 
dj = 2\ai\\ aj \. 



C. Distribution of Entanglement 

One way to study the time evolution of the spatial 
distribution of entanglement between the origin site and 
other sites in the spin chain is to look at the second order 
spatial moment of concurrence about the point of origin, 



Af(t) = ^|x 2 CW(t)| 



x^0 

where Co x (t) is the concurrence at time t, between the 
original site of the excitation, and the site at position x. 
C-Ox (t) can be replaced with the earlier expression for the 
concurrence when only one excitation is acted upon by a 
spin preserving Hamiltonian. 

M{t) = 5> 2 (2|aoW|K(t)|)| , (2) 

x^0 

= 2\a (t)\W(t) , (3) 

where 

W(f) = 5> a M*)l • (4) 

x^O 



In such a spin chain, with the additional constraint that 
couplings obey J*' l+1 = J y ' t+1 — T, for all sites i in the 
chain, then the overall z— component of spin is conserved. 
Here we will take r = 1, which will give general results 
up to a rescalling of t. For a spin chain of N sites, a 
basis can be formed from N vectors, |l)...|iV), with \k) 
corresponding to a chain with a single excitation at site 
k. Any Hamiltonian, if of the above form in this basis, 
corresponds to a symmetric tridiagonal matrix. Further, 
if the coupling parameters J\ — 0, then the main diagonal 
will also vanish, giving 



Hij = 8(i,j-l) + 8(i,j + l) 



H 



L + R , 



where Ljj = 5(i,j — 1), Rij — 5(i,j + 1). Using this 
expression for iJ, the evolution operator, U{t) can be 
found. 



U(t) = e 



_ p ~iHt _ e -i{L+R)t 



00 1 

j2 -(-itr(L+Rr , 

n=0 

00 n 



m 7~>n—m 



n=0 
00 



— ' (n — m)!m! 

m=0 v / 



L m R 



1 



[n — m)!m! 



2m 



where = R k , and = L. The only non-zero entries 
in Dk are along the k th sub-diagonal, which are all 1. The 
family of matrices, {-Dfc}!^, form a basis for U(t) and 
one can set, 



u(t) = c -MVx , 



1 



/=0 



{l + x)U\ 



= H)^£(-i)W- 



where 21 = n — x, 



1 



(l + x)W 



= (-i) x J x (2t) , 
where Jk is the k th order Bessel function of the first kind. 

1. Upper bound on W(t) 



D. Spin Dynamics: Ordered Chain 

A general Hamiltonian for a spin chain with near- 
est neighbour interactions and uniform couplings can be 
written as: 



J. 



In the case of the ordered chain we have not been able 
to obtain exact analytic expressions for but instead 
we can find bounds for W(t). Using and the approx- 
imation, 



Jx(z) 



for z > \x\ , 
'£cos(z-^-f) , for z < \x\ 



(5) 
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an upper bound for W(t) can be found, 

+00 +00 

w(t) = * 2 M*)I = E AM*t)\ . 

x — — 00 x— — 00 

+00 

= 53 \x 2 uzt)\ ■ 

X — — OO 

This can be simplified using the identity xJ x {a) 
1(^-1(0) + Ja;+i(o)) to obtain 

W(t) = \t 2 (Jx-2(a) + 2J x (a) + J x+2 (a)) 

X — — OO 

+t(J x - 1 (a) - J x+ i(a))\ , 
= (4t 2 + 2t) J2 \ J *( 2t )\ > 



1 



< (At 2 + 2t)\ — (At) 
V 7rt 

16 5 , 
» — =i 2 for t > 1 . 



S. Lower bound on W(t) 

A lower bound for W(t) can also be found by observing 
|J 2 fc(2t)| > J 2fe (2t) , 

OO OO 

(2k) 2 \J 2k {2t)\ > (M) 2 J2k(2t) , 

k— — 00 fe=— 00 

OC OO 

AJx(2t)\ > (2k) 2 \J 2k (2t)\ , 



X— — OO 



k— — c 



00 



> (2^) 2 J 2fe (2t) 

fc= — OO 



Using the identities Efel-oo( 2fc ) 2 J 2fe(a) = % |23|. and 
0, this can be rewritten, 

W(t) > 2t 2 . 

Since 2t 2 < W(t) < -7M for f > 1, then W[t) grows 
as f\ with 2 < A < 2.5. 



5. Approximation for M(t) 

Using the previous approxima- 

tion 101, which can be rewritten as, 

J„(z) « J^ z ( cos ( z - |) cos(^) + sin(z - f ) sin(f )), 



an expression for W(t) can be found. 
W{t) = £ x 2 \J x (2t)\, 



X— — OO 

2t 

E 



.' - \j — I cos(2t 

irt 2 



4>" 



By taking the average of | cos(2i — ^ — f )| over a time 
period from i = tf to i = tf + n, the time averaged value, 
~, is obtained. Using this as an approximation for the 
cosine term, the expression for W(t) reduces to: 



W(t) 



4 

7T V 7rf 



— 2t 

-Y- 



x=0 



4 [T2t{2t + l)(4t+ 1) 



7T V Trt 6 
32 5 

— > — 3-t 2 , as £ — ► 00 

37T2 

This can be used in the expression for M(t), 

M(t) = 2\a (t)\W(t) , 

n~2 32<l +24ti +4i2 

V 7Tt 7T 37rf 

128 2 

— > r— 3^ 7 a s t — > 00 . 

37T J 



E. Spin Dynamics: Disordered Chain 

It is well known that disorder, which in this case we 
take to be random J z coupling strengths, or random mag- 
netic fields in the z— direction at each site, leads to spatial 
localisation of the wave function |24j . In a spin chain with 
a disordered region centred on the initial site of excita- 
tion, the wave function will be partially localized, and so 
the amplitude at the initial site will fall off at a lower 
rate than in the case of an ordered chain. We now model 
the emission of an excitation from the centre site out 
through a disordered region and into (on either side), a 
semi-infinite ordered spin chain. 

To consider this we first note that the dynamics of 
an ordered semi-infinite spin chain, with an initial single 
excitation at the first site, can be solved in a similar way 
to the infinite spin chain. For an ordered semi-infinite 
spin chain we have 



U(t) 



-i(L+R)t 



00 1 

Y -(-unL+Ry 



(6) 



n=0 



For a semi-infinite chain, however, [L, R] ^ 0. 
(LR)i t j — 5i_j, but (RL)i_j = Si.j(l — If the initial 

excitation is at the start of the chain, then the resulting 
wave, after time t, will be the first column of U(t). Tak- 
ing the first column matrix elements of U(t) to be the 
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wave launched from site x = 0, as U(t) x ,o, and we can 
obtain, 



U(t) x , = J2 



{~it) n ( ( n 



i i i n—x / l n— x— 2 



= Hr(j K (2t)+ j x+2 (2*)) , 

resulting in a wave function with elements a x (t) (|T|). 
Ox(t) = (-i) x (J x (2t) + J x+2 (2t)), 

We now suppose that the random region, extending from 
x = —L to x = +£, emits an excitation at +L, into the 
ordered semi-infinite region with amplitude /(f). Shift- 
ing coordinates to this interface by setting x' = x— L, and 
then (for convenience) dropping the prime, we can exam- 
ine the wave propagation into the ordered semi-infinite 
chain (now x > 1) to be: 



ax +L (t) = J /(f)(-i)-i(__ 7 )J x (2(i-OK 



We now break the emission processes into two time bins. 
For t > h: 



a x+L (t) = J fV)(-ir-Hj—p)JxMt-*)W 
+ f t W)(-i) x -H^-p)M%t - t'Wt , 

= {-iY-\-^-)Jx{2{t-T)) f 1 f(t')dt> 

+ f mHT-^j^WHt - t'))dt' , 
= {-iY-\-^-)Jx{i{t-T))(3 

t — T 

+ f mHr- l ( T ^- F )M2{t - . 

(7) 

For large t\, we will neglect the second term in (J7J, 
which is equivalent to letting f(t) — (38{t — r). This 
corresponds to a simple model for the emission in which 
an excitation is emited, with amplitude /3, at time t = 
t <C t%. For t < r, there is no emission and a x +i(f) = 0, 
while for t > r: 

a x+L (t) = /3(- l r- 1 (-^)J :c (2(i-r)) . (8) 
t — r 

The second order moment of concurrence can be split 
into two parts, M {t) and Md(t). M (t) is the contribu- 
tion to M(t) from the ordered region, and M<z(t) is the 
contribution from the disordered region. For t > t: 



M, 



(t) = l\ aQ {t)\Y,{x + LY^\PJ x {2{t-T))\, (9) 



= 2-\p—(\cos(2T--)\Y(2x + L) 2 2x 

V X — 1 

T 

+|sin(2T- ^)\"£(2x + L-l) 2 (2x-l)) , 

x=l 

(10) 

where T — t — r. 

Since the average value of | sin(z)| = | cos(z)| = -, the 
above can be approximated by 

Wo{t) = ^rT(x + L) 2 X, 



X — l 



,32L 



^-V(8T 4 + (— - 8)T 6 + {AL 2 -8L + 3)T 2 

Tz 7T2 O 

+(4L - 2 - 2L 2 )T - (18 + ^— + 2L 2 )) , 
161, 



-(t-r)^ , fort>0. 



From © , to obtain M a (t) , we require the amplitude of 
the wave at the origin or initial site. Since the probability 
of the excitation entering the semi-infinite chain to the 
right and left is |/3| 2 , the probability of the excitation 
remaining within the disordered region is 1 — 2\j3\ 2 . If the 
probability at the initial site is proportional to this, then 
the amplitude of the wavefunction at the initial origin 
will be 7v /l - 2|/3| 2 , with | 7 | < 1. 



1 27T2 6 

+ (AL 2 -8L + 3)T 2 + (4L - 2 - 2L 2 )T 
-(18+^ + 2L 2 )), 



32|/3 7 yi~2l 

3 

7T2 



-(i-r)3, fori»0. 



(11) 



There is only a finite disordered range, so Wd{t) is 
bounded from above by 2L 2 . This in turn leads to an 
upper limit on M d (t) of 4| 7 yl - 2\/3\ 2 \L 2 . Since M Q cx 
t§ and M d < 2/3L 2 , 



x=l 



M(t) « M (t) « - r)l, far * > . 

7T2 

(12) 

In summary, we have examined the emission of a sin- 
gle excitation from a central site on a spin chain in the 
two cases where the chain is ordered and where the ori- 
gin site is surrounded by a finite region of disorder. We 
have found that the spatial expansion of the two-site en- 
tanglement, between the origin site and other sites on 
the chain, is much more rapid in the latter case where 
disorder is present, than in the ordered case. Indeed the 
second order moment of the spatial extent of the con- 
currence expands super-ballistically. This effect could 
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be used to rapidly distribute entanglement throughout 
regions of a spin chain. If, rather than distribute en- 
tanglement throughout a chain, one wished to spatially 
separate an entangled state one might instead appeal to 



a recent scheme outlined in 14]. 
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